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1. Introduction 


The Global Ozone Monttoring Experiment (GOME) was launched on ERS-2 i April 1995 and 
has been continuously operational since then. It is a scanning nadir-viewing spectrometer, with its 
primary scientific objective being to retrieve total column ozone globally. A more detailed de- 
scription of the instrument can be found m /17/. In common with previous instruments to measure 
total column ozone from space such as TOMS and SBUV, it measures the back-scattered radiance 
from the Earth's atmosphere and surface, and the solar irradiance which is viewed via a difluser 
plate to provide a reference spectrum at comparable intensity. These difluser plates have been 
found to be subject to degradation (see /2/ for example) particularly when subject to shorter 
wavelength ultra-violet heht, and efforts have been made to characterize this degradation for in- 
struments such as SBUV/2, where the diffuser plate was exposed for a total ofaround 750 hours 
between 1979 and 1986. GOME has been designed with a cover for its diffuser plate in an attempt 
to minimise this degradation, with exposure usually being for a short time for one orbit each day 

to obtam a reference solar spectrum, and characterisation of any degradation is possible by means 
ofthe on-board PCr/Ne calibration lamp. This technical note investigates to what extent, if any, 
GOMEYs difluser plate has degraded, and to see ifthe measures taken have reduced the eflect. 


The time period for the analysis was a little more than one and a half years, using the monthly cal- 
ibration data from June 1995 until January 1997. During the course of this analysis, a problem 
was discovered with the algorithm used for the calculation of the difluser reflectivity, which 

meant that trends in the dark signal were causing difficulties i the interpretation of the reflectivity 
data. Consequently the two components ofthe dark signal were also analysed, and the results not 
only allowed a better understanding of the diffuser behaviour, but also provided valuable informa- 
tion on the detector degradation over this same period of time. 


2. ESRIN QA Tools and Data Set 


The mann tool used for this investigation was the Extended Rascals for GOME (ERGO) system 
produced under ESA contract by Dornier and SRON (/37, 47). This is a software suite designed 
to provide quality assurance measures about GOME on both an automatic and interactive basis. 


Every month a set of calibration files are generated which contam extra information about the in- 
strument, as obtamed from the monthly calibration sequence. Amongst these are measurements of 
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difluser reflectivity calculated once per spectral line per channel per orbit, as obtained using an al- 
gorithm specified in /5/ and updated in /6/, and the two dark signal components. These algo- 
rithms are described further below in $3. 


The data set consists of some 20 monthly calibration sequences which are not quite regularly 
spaced in time - data were unavailable for some months as well as there bemg a couple ofocca- 
sions on Which the calibration sequence was performed more frequently. The table in Appendix A 
gives information about the orbits available for analysis from each monthly calibration sequence. 


3. Algorithm Descriptions 


The algorithms used for the calculation of difluser reflectivity and dark signal components are 
summarised in the following subsections. 
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3.1 Diffuser Reflectivity Algorithm 

The algorithm used by ERGO to calculate difluser reflectivity may be broken down into the fo- 
lowing components: 

e Find all difluser products for an orbit 

e Find all lamp products after the last diffuser product 

e Tonore the first lamp products where the integration times are all minimum 

e Locate the centre ofeach line 


e Sum five pixels around this line centre 


e Calculate baseline from average of the minima of ten pixels either side ofthe centre 
e Subtract this baselne 


e Normalise according to integration time 

„ Average the difluser measurements 

„ Average the lamp measurements 

e Calculate the lne reflectivity as the ratio ofthese averages. 

The wavelengths and corresponding pixel numbers ofthe calibration lamp lines used are shown in 
the tables in Appendix B. As noted previously, the algorithm as it stands yields one ratio per line 
per orbit, and no error information on this measurement. Some account of lamp power stability in 


taken by not using measurements for the first 130s after it has been switched on, thereby avoiding 
its power-up phase. 


3.2 Dark Signal Algorithm 


e Find all dark products for an orbit 
e Get their corresponding integration times 
e Lmear fit the integration times against the dark signal, notmg Band 1A integration status 


e Return the ofŝset as the fixed pattern readout noise and the gradient as the leakage current 


4. Preliminary Analysis of Diffuser Reflectivity Data 


A preliminary analysis was performed on the ERGO data in July to see what conclusions could be 

drawn from the data set as it stood at that time. It can immediately be seen (Figure 1, topmost 

plot) that the data come in clusters of several points, corresponding to each monthly calibration 

sequence. This period is not completely regular since there were occasions (March and June 1996 

for example) on which more than one monthly calibration occurred, and occasionally the se- 

quence was not performed as scheduled on the 28 8 ofeach month. The calibration is performed 
over the Pacific Ocean, so only the first two of the calibration orbits at most could be aflected by 

the South Atlantic Anomaly (SAA), an area of higher than normal background radiation, and then 

only during the daytime when the orbit track passes through this area. 


There are two measures Which may be taken in an attempt to identify statistically erroneous points 
Which can then be removed, yielding a smoother data set for further analysis. The first of these is 
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to obtain the mean reflectivity ofall the poimts and discard those lying more than 3 standard devi- 
ations (SD) from the mean (see Table 1 below); the second is to analyse the spread of the points in 
the individual clusters and to discard those clusters where this spread is significantly (more than 

2.5 SD - since we are dealing with a group of lines rather than individuals the selection criterion is 
stricter) greater than the others (see Table 2 below). The result of applying these transformations 

to the first lne mm channel 1 can be seen in Figure 1, where the points to be discarded are marked 
with a star. 


Channel Nov95 Early Feb 96 Mar 96 Apr 96 May 96 Jun 96 
One - - - - 1-7,9, 3 

Two - - - 2-6,8 1,7 - 
Three 6 9 - 1-4,8, 11, 13-18 5-7, 10, 12 10 
Four 1 - 20, 21 2,4-6, 8-21 1-4,7 1-3 


Table 1: Bad Line Occurrences Listed By Line Number In Channel 


Channel Dec 95 mid Mar 96 end Mar 96 Apr 96 May 96 Jun 96 
One - 1,8 - - 9 5-7 


Two - 6 - 7 2,4 1,3,8 


Three - - - - 1,9, 11, 13 4-7, 12 
Four 21 - 16-19 23,1] 8; 12 - 


Table 2: Bad Cluster Occurrences Listed By Line Number In Channel 


From these two tables it can be seen that there is no consistent pattern where there are problems in 
all channels for a particular month, except for May 1996. There are some 100 occurences listed 
here, ofwhich 40 are related to one particular orbit (5346) m April 1996. Although this orbit does 
pass close to the SAA, it would be unreasonable to assume that this was the cause ofthe bad data, 
since other orbits which are similarly close do not show similar behaviour. For example, orbit 

4488 which is the first orbit ofthe calibration sequence in late February 1996 has no bad data 
points associated with it. Further to this, ofthe remaining 60 bad pomts only 26 occur in orbits 
where they might have been affected by the SAA. This result is also expected since the long inte- 
gration times used for the difluser characterisation, which would be sensitive to perturbations in 

the background radiation environment due to the low signal, are time-lned to occur durmg the 
nighttime part of the orbit. 


Once all the data had been processed in this way, it was possible to examine dependencies of the 
difluser reflectivity on wavelength (by lme and by channel) and time (orbit). A surface plot for 
channel 1 is shown in Figure 2, where dependencies on both waveleneth (higher reflectivity at 
shorter waveleneths) and time (reflectivity mcreasing with time) can be seen. 
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When an average is taken over wavelength for each channel, only the time dependency remajns, 
and a straight lne fit to the data gives the following table, where all reflectivity data have been 
multiphed by one thousand. The data are shown in Figure 3. 


Channel No. 1 2 3 4 
Mean 830 570 650 540 


Yo | Year 14 -0.9 1.8 1.1 


Table 3: Preliminary Diffuser Reflectivity Values “ 1 000 000 


The conclusions one can draw from this are rather limited, since the data are obviously still noisy, 
but essentially there is no discernible trend in channels 2, 3 and 4. Channel 1 however shows an 
mncrease in reflectivity of around 1499 a year, with the reflectivity bemg higher at shorter wave- 
lengths (as seen in Figure 2). This was rather unexpected behaviour (if anything a decrease in re- 
flectivity at short wavelengths was predicted /2/) and so eflorts were made to a) reduce the 
noismess ofthe data, particularly for channel 1 and b) find an explanation for the observed behav- 
lour. 


5. Diffuser Reflectivity Algorithm Enhancement 


Given the large number of direct lamp measurements per orbit, and the fact that several measure- 
ments via the difluser are also made, combining these after averaging into just one ratio with no 
error information as was done previously is obviouslY resultmg in some loss of information. Con- 
sequently it was decided to ratio all lamp measurements agamst all lamp via diffuser measure- 
ments, combining these to provide a mean and variance. This error measure then allows the 
measurements for individual lines to be combined in the standard way /7/ which gives more 
weight to those with a small error i.e: 
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thus giving an optimal estimate of channel reflectivity. The effect this technique has can been seen 
by examining Figure 4, which shows the errors for each cluster for the individual lines of channel 
1, using a different symbol for each cluster. Correlation is obviously low between error size and 
cluster (im each cluster one finds lines with low and high errors, and this pattern is different for 
each cluster), whereas 1t is noticeable that lines 1, 2, 3 and 6 have consistently higher errors. This 
means that when the data are combined to provide a channel reflectivity, these lines will provide 
only a very small contribution. 


Figure 5 shows the behaviour with time of the reflectivity calculated for individual lines from 
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these values will decrease for the channel average using the above method. 
The optimal estimation technique was included in an upgrade (August 1996) to the ERGO sys- 


tem, and all the monthly calibration data were re-processed, together with later data from July, 
August and September 1996. They were again filtered statistically as described above in $4 (note 
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that the revised algorithm makes little difference to the individual values, so the same lines were 
removed, but there is now an error estimate for each line). The corresponding table to that shown 
i the prelimmary analysis section is shown below, and the data are plotted as the bottom graph in 
Figures 6 to 9. 


Channel 1 2 3 4 
Mean 620 560 540 540 
Vo | Year 92 1.4 -1I7 -0.9 


Table 4: Revised Diffuser Reflectivity Results “ 1000 000 
Although both the noise and the size of the trend m the data for channel 1 has been reduced and 


brought more into line with other results, it is still significantly different and so efforts were 
turned towards finding an explanation for this difference 


6. Dark Signal Analysis 


During an early review of this document, C. Caspar (ESTEC) noted that there was a problem with 
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tribution to the total beg underestimated, and the line signal overestimated. Under circumstanc- 
es where the dark signal is much greater than the line signal, such as occurs when a very weak 
lamp line is viewed via the difluser, it is possible for the difluser reflectivity to be substantially 
overestimated. Consequently, an analysis ofthe dark signal was performed which allowed the ef- 
fect to be quantified. 


The dark signal for GOME is defined as bemg comprised oftwo parts - a constant value of be- 
tween 140 and 150 binary units (BU) which is the fixed pattern readout noise (FPRN) and a time 
dependent component of around 2 binary units per second which is the leakage current (LC). 


Results for January 1997 for two ofthe four channels are shown in Figure 10, using only those 
products where the integration is complete for all bands so as to avoid a known cross-talk problem. 
Immediately apparent im the FPRN is a superimposed spectral signature probably due to light 
reflected from a baffle when the muror is pointing to deep space. 


To investigate this further, the data were decomposed into the components from each different 
integration time. An example of'this is shown from June 1996 for channel 4 m Figure 11. Looking 
at all the data, it became clear that this spectral signature is only present for integration times of35s 
and 24s, where the monthly calibration sequence performs dark signal measurements in the 
daytime part of the orbit. 


A further problem was also seen with the measurements in channel 2 at 1.5s, channel 3 at 0.375s 
and channel 4 at 0.09375s where a slow change i the dark signal across the detector elements 1s 
seen. The problem of some dark signal readouts beg dependent on the integration status of other 
channeles was noted in the functional performance tests performed at the TNO Institute of Applied 
Physics (TPD) in Holland on the breadboard model (BBM) /8/. However, the monthly calibration 
timeline does not vary, and so it was sufficient to identify where these eflects occurred. They are 
shown for the example case of June 1996 i Figure 12. Consequently measurements at these 
problematic integration times were excluded from further trend analysis. Note that discardig the 
daytime measurements removes any potential problems from South Atlantic anomaly eflects, since 
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The analysis ofthe remaming data was performed by taking the middle file ofeach monthly 
calibration sequence, calculating the gradient (LC) first from all valid integration times, and then 
using this to extrapolate from the shortest integration time to get the ofiset (FPRN). Trends were 
calculated for both the FPRN and the LC, and for the noise on these measurements. The results 
are shown in Table 1 below. 








Ch. FPRN Notse LC Noise 
1 H0.6 Sub, 13,3 434 
2 =l,2 H14.9 H14.9 435 
3 0.3 8.6 [33 14733 
4 H0.6 E, 17.8 88 














Table 5: GOME Dark Signal Trends; Yo per year 


From this one can conclude that there is no significant change in the fixed pattern readout noise, 
but that a lnear increase in the leakage current is seen, Which is also becoming noisier with time. 


7. Impact of Dark Signal Trends on Diffuser Reflectivity Analysis 


To calculate the impact of the dark signal trends, the dark signal measurements were analysed at 
the positions of all the lamp lines used for diffuser characterisation, and the values ofthe mean - 
minimum difference calculated for each monthly calibration sequence. These values were then 
used to adjust the lamp measurements via the difluser, and the reflectivity values recakculated, ne- 
glecting data from lines so weak that they were less than the dark signal -- 3SD. In the table below 
are shown the data for the lines m channel 1, and it can be seen immediately that the effect of the 
dark signal trend is to cause both the absolute values of all lnes viewed via the diffuser, and the 
trends therem, to be overestimated. 


Approx. Annual Apparent Start / End Intensity Intensity 
Line Intensity Change Intensity Mean - Mn Overest. Overest. 
(BU) Mean - Mn Change Start End 
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1 40 10 2599 23 36 
2 85 16 1999 19 46 
3 Lo ko 2099 24 49 
4 190 16 ZI 20 40 
5 120 14 1299 19 38 
6 2 3 So 22 21 
1 205 9 490 19 35 
Ŝ 230 14 690 18 41 
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Approx. Annual Apparent Start / End 
Line Intensity Change Intensity Mean - Mn 
(BU) Mean - Mn Change 
9 319 8 2Y9 29 44 
10 1000 20 2Yo 15 49 
l 445 9 2Yo 23 32 


These corrections were performed to all the data used, and the difluser reflectivity values recalcu- 
lated. The results and the zero change lines are shown i Figure 13. On these plots are also shown 
the errors associated with each point (3 SD, and quite large since they involve the combinmg of 
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8. Comparison with Other Results 


TPD and ESTEC presented therr results at the GOME QA workshop held at ESRIN (Frascati, Ita- 
ly) in early September 1996. Therr values, reproduced with permission, are given in the table be- 
low and, together with the ESRIN results after dark signal correction all given to two significant 
figures. The margins given on the 90 / year figure are 1 SD in the straight line fitting coefficients 

and are as plotted im Figures 6 to 9, except for the ESRIN data which are those as shown in Figure 
13 (after dark signal correction). Note that i these plots the x axis is in days since launch for TPD 
and ESTEC, and orbits since launch for ESRIN. 


Ch Mean SD 00 / Year 
TPD - ESRIN ESTEC TPD ESRIN ESTEC TPD = ESRIN ESTEC 
1 570 570 590 30 14 12 25116 -3.011.0 -2.0-10.5 
2» 520 520 500 20 14 9 19-49 1.040.5 2.5:40.5 
3 470 520 480 30 9 Di 6E16  0.012.0 0.5-E0.5 
4 540 530 540 20 4 S - 168 0.010.5 0.5:H0.5 


Table 6: Comparison of Diffuser Reflectivity Results ~ 1000 000 


Comparing these results it can be seen that there are the greatest margins on the TPD data set 
Which is due partly to the fact that it contams fewer (6) samples. The larger margins on the ESRIN 
data set are for the reasons given at the end of section 7. It can be seen that agreement is generally 
good, particularky between the ESRIN and ESTEC data sets which cover a much longer period 
than that of TPD. 
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9. Summary 


GOME difluser reflectivity data as calculated using the ERGO algorithm are iherently noisy. Op- 
timal estimation as described above provides a good way to combine them, making full use of the 
data available at ESRIN. Channels 2, 3 and 4 show insignificant changes over the 80004 orbits for 
Which they have been monitored at ESRIN, and agreement with results from TPD and ESTEC 1s 
good for these channels. 


Channel 1 showed a significant increase in difluser reflectivity, but this was shown to be due to a 
weakness in the ERGO algorithm. It affected Channel 1 predominantly since there are many weak 
lamp lines used in this region for difluser characterisation. 


Trends in the dark signal were mvestigated and shown to be capable of explaining all apparent 
trends i difluser reflectivity. 


10. Conclusions and Recommendations 


Over a period of more than eighteen months, the following conclusions regarding trends in the 
GOME dark signal and difluser reflectivity have been reached. 


e No sienificant change is seen in the fixed pattern readout noise 

e There is an increase of~1590 per year in leakage current for all detectors 
e Leakage current measurements are becoming much noisier with time 

e No significant change is seen in the difluser reflectivity in any channel 


These results, from in- flight measurements, can be utilised to predict the likely degradation which 
will be encountered in missions with similar components, such aa SCIAMACHY. One can further 


sencludn fate MEANS SAKENdP LIRAMLMIbs Asreatr dom degradation have been 


It is recommended that the ERGO system calculate the dark signal separately from the monthly 
calibration files, and use this for a proper dark signal correction of the diffuser. Ifthis proves 
impossible, significant improvement would still be gamed by using only the strongest lines for the 
calculation, since this is where the impact of the dark signal trend is least. An alternative method 
would be to include a parameterisation ofthe dark signal behaviour based on the results of'this 
analysis, and to use this to correct the results produced by the existmg algorithm, though this 
approach would be the least satisfactory. It is not possible to make a satisfactory dark signal 
correction based on measurements made whilst the lamp is on. 
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Figure 1. Statistical Filtering of Data from First Line in Channel 1 
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Figure 2. Surface Plot of Channel 1 Reflectivity Data 
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Figure 3. Prelimmary Results for Mean Channel Reflectivity 
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Figure 4. Errors from Each Cluster for Each Line of Channel 1 
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Figure 5. Behaviour over Time for All Lnes of Chamnel 1 
(the numerals correspond to those in Appendix B) 
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Figure 6. Channel 1 Reflectivity Comparison - No Dark Signal Correction for ESRIN Data 
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Figure 7. Channel 2 Reflectivity Comparison - No Dark Signal Correction for ESRIN Data 
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Figure 8. Channel 3 Reflectivity Comparison - No Dark Signal Correction for ESRIN Data 
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Figure 9. Channel 4 Reflectivity Comparison - No Dark Signal Correction for ESRIN Data 
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Figure 10: ERGO Dark Signal and Solar Spectrum, January 1997 


Figure 11: Dark Signal Measurements for Channel 4, June 1996 
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Figure 12: Integration Times and Channels where Variation of Dark Signal Across Detector 
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Elements is Seen 


Figure 13: Diffuser Reflectivity (x 10,000) for all Channel, Dark Signal Correction Applied 
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Appendix A 


Monthly Calibration Data Sets 


— kM UA 
1 27 June 1995 67 965 - 967; 969 (4) 
2 28 July 1995 98 1410 - 1413 (4) 
3 28 August 1995 129 1854 - 1857 (4) 
4 28 September 1995 160 2298 - 2301 (4) 
5 28 October 1995 190 2726 - 2730 (5) 
6 28 November 1995 221 3171 -3174 (4) 
ĝi 28 December 1995 251 3600 - 3604 (5) 
8 04 February 1996 289 4144 - 4148 (5) 
9 28 Februrary 1996 313 4488 - 4491 (4) 
10 13 March 1996 327 4684; 4687 (2) 
11 26 March 1996 340 4874 - 4878 (5) 
12 28 March 1996 342 4902; 4904 - 4906 (4) 
13 28 April 1996 3173 5347 - 5350 (4) 
14 28 May 1996 403 STT6 - S780 (5) 
15 28 June 1996 434 6220; 6221; 6223 (3) 
16 28 July 1996 464 6649; 6650 (2) 
17 28 August 1996 495 7092 - 7096 (5) 


18 28 September 1996 526 7536 - 7540 (5) 
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Lime 
Number 


19 06 November 1996 565 8094 - 8098 (5) 
20 28 January 1997 648 9282 - 9286 (5) 
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Appendix B 


Lamp Lines Used For Diffuser Calibration 


Channel 1 Channel 2 
Wavelength Pixel Lme Wavelength Pixel 
/nm Number Number /nm Number 
244.08 313.79 1 321.91 215 MI 


248.19 353.41 2 3247 368.00 


2 262.88 415.23 3 337.92 415.96 


4 266.02 503.16 4 39241 542.43 
5 273.48 569.93 o 369.53 698.64 
6 281.03 638.76 6 312.82 728.42 
1, 28341 657.69 di 390.99 893.48 
8 293.06 749.28 8 392.03 903.05 
9 299.88 812.66 
10 304.35 854.03 
11 306.56 874.64 
Channel 3 Channel 4 

Line Wavelength Pixel Lme Wavelength Pixel 
Number /nm Number Number /nm Number 
1 425.55 145.6 1 588.35 44.42 
2 427.60 1552 2 594.65 122 
3 429.09 162.2 3 597.72 85.92 
4 437.25 200.6 4 603.17 109.43 
2 460.20 309.8 5 607.60 130.53 
6 492.36 464.5 6 609.79 140.45 
7 503.92 520.3 7 613.01 15521 
8 540.21 694.7 8 616.53 171.30 
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Channel 3 Channel 4 
Line Wavelength Pixel Lme Wavelength Pixel 
Number /nm Number Number /nm Number 
9 556.43 772.0 9 621.90 196.02 
10 574.99 859.4 10 626.82 218.84 
Hl 576.60 866.8 1 630.65 236.66 
12 580.61 885.5 12 638.47 21323 
13 582.18 892.7 13 653.47 344.13 
14 588.35 921.3 14 660.08 375.64 
15 594.65 950.4 15 668.01 413.61 
16 so. 2 964.1 16 693.14 534.83 
17 603.17 988.6 l i 717.59 653.62 
18 607.60 1008.5 18 724.72 688.28 
19 744.09 782.40 
20 749.09 806.62 


sal kodo 829.28 
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